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Алгоритм моделирования искусственной 
иммунной системь с оператором селекции Саати 
и одномерньм локальньм поиском 


Розглядається метод умовної оптимізації у безперервному просторі багатьох змінних, заснований 
на моделюванні імунної системи людини. Запропонований новий оператор селекції клітин для клонування, 
що базується на методі аналізу ієрархій (МАТ) Сааті. Для інтенсифікації роботи евристичного алгоритму 
використаний оператор одновимірного локального пошуку. 
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Рассматриваєтся метод условной оптимизации в непрерьтвном пространстве многих переменньтх, 
основанньшй на моделировании иммунной системь человека. Предложен новьгй оператор селекции клеток 
для клонирования, основанньй на методе анализа иерархий (МАЙ) Саати. Для интенсификации работьт 
звристического алгоритма использован оператор одномерного локального поиска. 

Ключевьге слова: иммунньге системьк, оптимизация, Саати, локальньгй поиск. 


Вступ 


Чимало інженерних, соціальних та економічних задач можуть бути сформульо - 
вані у вигляді задач оптимізації деяких цільових функцій багатьох змінних. При цьому в 
більшості випадків цільові функції виявляються мультимодальними. Наявність кількох 


І59М 1561-5359 «Штучний інтелект» | 2013 Мо 4 101 


Желдак Т.А., Слєсарєв В.В. 
ЮНГ б б 


локальних оптимумів цільової функції одночасно з одним чи декількома глобальни- 
ми оптимумами традиційно розглядається як суттєвий недолік. Існує цілий ряд мето- 
дів 1 алгоритмічних прийомів, спрямованих на вихід з локальних оптимумів з метою 
досягнення єдиного глобального. 

Задача ускладнюється, коли йдеться не про один, а про декілька критеріїв та мно- 
жину рішень, недомінованих одне одним, серед яких все одно намагаються віднайти 
єдиний оптимум. В реальних застосуваннях часто існує необхідність пошуку не єдиного 
оптимального розв'язку задачі, а їх сімейства. Кожне з отриманих рішень, що в загально- 
му випадку є субоптимальним, дозволяє розглянути одразу кілька можливих результатів. 

Таке багатоваріантне розв'язання задачі дозволяє отримати метаевристика, відома 
як метод моделювання штучної імунної системи. 

Запропонований нижче імунний алгоритм (ТА) імітує властивості природної імун- 
ної системи і заснований на принципах соматичної теорії (1 | і мережевої гіпотези |2|. 
Соматична теорія стверджує, що збільшення різноманітності антитіл відбувається за ра- 
хунок соматичної рекомбінації і мутації генів. В рамках мережевої гіпотези обгрунто- 
вується припущення, згідно з яким контроль розмноження клонів здійснюється в резуль- 
таті взаємного розпізнавання антитіл, що функціонують як єдина мережа | 11, (31, Г41. 

В процесі еволюції імунна система вищих ссавців набула здатності знищувати 
антигени за допомогою антитіл (рис. 1). 
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Рисунок 1 - Механізм формування антитіл в імунній системі 


Клітинна 
пам'ять 


На рис. 1 показана робота імунної системи з точки зору процесів обробки і управ- 
ління інформацією. Користуючись цією схемою, інформаційні процеси в імунній системі 
можна представити у вигляді алгоритму, який вкладається у наступну послідовність 
формалізованих кроків, докладно розглянутих у |З: 

1 Розпізнавання антигену. Цей процес відповідає встановленню вигляду задачі 


оптимізації / « / (х 7- тіпо / ра ), де п- розмірність вектора Х дійсних змінних 
Херев" 


х, Є Х , заданих на області припустимих значень Р. 


2 Генерація антитіл клітинами пам'яті. Це відповідає пригадуванню успішного 
вирішення аналогічної (або схожої) задачі в минулому. 
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3 Обчислення афінності (корисності) клітини Х,, що є еквівалентним пошуку 
оптимального рішення: ф, 2 -/ (х 4) для задачі мінімізації. 


4 Диференціювання лімфоцитів: збереження відповідного рішення для наступного 
кроку пошуку та видалення зайвих кандидатів на рішення. 

5 Стимулювання 1 придушення антитіл в залежності від афінності. Дослідження 
локальних оптимумів супроводжується підтриманням розмаїття напрямів пошуку. 

6 Розмноження антитіл з використанням пам'яті та випадкових генів. 

У природі описаний процес є безперервним - антитіла продукуються організмом 
як до відомого антигену, так і до його можливих мутацій. При моделюванні процесу 
на комп'ютері роботу алгоритму можна зупиняти після певної кількості поколінь або 
ж через заданий наперед час. 

Результатом роботи алгоритму буде сімейство рішень, що відповідають антитілам 
у клітинах пам'яті, максимально пристосованих до антигену. У випадку пошуку екстре- 
мумів багатомодальних функцій це сімейство відповідатиме масиву з векторів-рішень, 
афінність яких максимально наближена до глобального оптимуму, тоді як самі рішення 
попарно відмінні не менш як на задану порогову величину. 

Показниками ефективності роботи алгоритму можуть бути час виконання та кіль- 
кість звернень до цільової функції при відшуканні глобального оптимуму, а також надій- 
ність (повторюваність) отриманих рішень при повторних запусках. Останнє є доволі 
важливим, адже евристичний метод, що розглядається, відноситься до класу методів 
випадкового пошуку і має збіжність лише за ймовірністю. 

Метою даної роботи є розробка такого алгоритму, що реалізовував би метод 
моделювання штучної імунної системи для вирішення задачі багатовимірної умовної 
оптимізації багатомодальних безперервних функцій та забезпечував при цьому макси- 
мально можливі показники ефективності. 


Аналіз існуючих алгоритмічних рішень 


Історично першою реалізацією ідеї моделювання штучних імунних систем була 
ідея селекції та витіснення незалежних пошукових агентів, запропонована у алгорит- 
мі СТОМАТ.С (5). Даний алгоритм, що є поєднанням головних ідей генетичних алго- 
ритмів та еволюційних стратегій, використовує двійкове кодування дійсного простору 
пошуку, селекцію певної частини (близько половини) кращих особин 1 їх рівномірне 
клонування, мутацію зворотно, пропорційну корисності, виключення найменш корис- 
них антитіл та стиснення популяції. 

Оператор стиснення популяції вилучає надлишковість популяції, виключаючи 
гірше з близьких рішень, що задовольняють вимозі 


Її Ха - Хь| «т» ас, абеД:М, |, (1) 


де г - певний поріг близькості; и - символ норми, яка розраховується в залеж- 


ності від метрики простору; М, - розмір популяції. 

До недоліків алгоритму С.ОМАІТ С слід віднести відсутність адаптації, обмежену 
точність, обумовлену двійковим кодуванням, та невикористання як такого механізму 
пам'яті. 

Інший алгоритм моделювання імунних систем - алгоритм ВСА (|6| - фактично 
є версією еволюційної стратегії з двійковим кодуванням і оператором стиснення. Він 
використовує одразу два оператори мутації: у 3/4 випадків - бітову, а у 1/4 випадків -- 
суміжну (кілька бітів, що стоять поруч, інвертуються одночасно). Метод показує високу 
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ефективність пошуку при невеликих розмірах популяції М, -3--5, однак, як і СТОМАІС, 


на багатовимірних тестових задачах за якістю рішень програє відомим еволюційним 
стратегіям. 

Певним кроком вперед є алгоритм орі-АїМеї, головною особливістю якого є робо- 
та безпосередньо у дійсному просторі та збереження всіх знайдених локальних 1 глобаль- 
них оптимумів цільової функції у пам'яті, розмір якої динамічно змінюється. Алгоритм 
використовує також динамічний розмір популяції, керований механізмом стискання |3|. 
Здійснюється мутація кожного клону, пропорційно корисності його батьківської кліти- 
ни за формулою 


ой здо сь во 
хо зо нау МОЇ, аеро" 14ЛМ, 0) 
де ХУ - і-й клон /-ї клітини попередньої популяції; с, - радіус мутації для 


клонів /-ї клітини; Ф - відносна афінність / -ї клітини у поколінні фіз их 8 р -ві- 
Фах " Фипіп 

льний позитивний параметр методу; М (031) - нормально розподілена випадкова вели- 

чина з математичним очікуванням 0 і середньоквадратичним відхиленням (СКВ), що 

дорівнює І. 

Практична реалізація даного алгоритму на ряді функцій |7| показала його знач- 
ну перевагу над алгоритмом СТОМАІТ С у швидкості отримання порівняного за якістю 
рішення. 

В роботі |4| запропоновано доповнити алгоритм орі-АїМеї оператором випадко- 
вого локального пошуку. Останнє робить деякі реалізації алгоритму порівняними за 
ефективністю роботи з відомими генетичними алгоритмами та еволюційними стратегіями 
на класі розв'язуваних задач. Не випадково - 1 у цьому головна цінність орі- ДїМеї -- 
він є основою для кількох більш сучасних алгоритмів, що реалізують викладений вище 
метод моделювання імунних систем. 

При вирішенні багатьох прикладних задач з використанням методу моделювання 
імунних систем застосовується алгоритм НІА, запропонований в (1 | 1 удосконалений 
в 7). Його можна розглядати як модифікацію методу орі-АїМеї. На відміну від останньо- 
го, розмір популяції в алгоритмі НІА залишається постійним в процесі всього ітерацій- 
ного процесу й підтримується оператором стискання (1). 

Крім того, кожна з клітин має додатковий атрибут /, званий віком клітини. За 
його допомогою визначаються так звані «елітні» клітини, що прожили / 21... 


Це - рішення, які тривалий час не були покращенні мутацією, тому вважаються 
локальними оптимумами. Вони заносяться в пам'ять й використовуються для генера- 
ції клонів у кожному новому поколінні. 

Мутацію клонів здійснюють двома способами. Перший, відповідальний за інтенси- 
фікацію пошуку в околі знайденого рішення, полягає у покоординатній зміні клона з ви- 
користанням адаптивного нормального розподілу рівня мутації 

Х9 -Х9 Ма), З) 
де а, (1) - СКВ нормального розподілу, що для кожної ітерації / визначається за- 
лежністю 
хе зх бе Ле )) є ГО) 
а, ( 1) пе / / / / й 
а я 1) оуегугі5е 


(4) 
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Другий спосіб має на меті диверсифікувати напрямок пошуку за рахунок обсте- 
ження ще не розглянутих ділянок 1 полягає у зміні однієї координати і випадковим 
числом з усієї області пошуку 


Хі - ОДБ Рол і зе Й : п| Ю (5) 
Туг п - розмірність простору; / рань сли ,) - рівномірно розподілена випадкова 


он, у. 


і, тах 


величина в межах Др з Їх ге; 


і,тіп 

Алгоритм НІА, у порівнянні з попередніми, демонструє надійність пошуку гло- 
бального оптимуму багатоекстремальних функцій за доволі обмежений час. 

Деякі дослідники, зокрема |8|, 19), включають до алгоритмів моделювання імун- 
них систем оператори локального пошуку, які базуються на наведених вище методах 
безумовної оптимізації. Втім, застосування методів сходження на гору, Нелдера-Міда 
та множників Лагранжа, зокрема, у методі 5ІА |8| для розв'язання задач з обмежен- 
нями, хоч Й посилює локальну збіжність алгоритмів, втім, не істотно впливає на ймовір- 
ність отримання глобального оптимуму. 

У алгоритмі Рорі-АїМе! |9| застосовано оператор копіювання координат, який 
працює наступним чином: якщо після мутації по якійсь із координат афінність клону 
покращилася, нове значення за цією координатою намагаються застосувати до решти 
вимірів. Такий оператор показує високу ефективність в разі симетрії функції по коорди- 
натах, властивої для більшості тестових функцій, наприклад, Растрігіна, Еклі, Швефеля, 
Гріванка та інших |2). Однак, при розв'язанні реальних задач, таких як навчання нейрон- 
них мереж прямого розповсюдження, чи основаних на радіально-базисних функціях, 
цей оператор абсолютно не ефективний. 

Більш доцільним бачиться введення у загальний алгоритм роботи штучної імунної 
системи оператора інформаційного обміну між пошуковими агентами (рекомбінації). 
Останній, відомий нам з генетичного алгоритму під назвою кросоверу, дозволяє значно 
покращити якість рішень 1 досягати потрібної якості в реальному масштабі часу за ра- 
хунок інформаційного обміну. 

Зокрема, у (10) пропонується для кожного з утворених клонів з невеликою ймовір- 
ністю Р. 20,2 виконувати його кросовер з випадковою клітиною пам'яті. При цьому 


СГО55 7 
для кожного з генів рівноймовірно, залишаться вони без змін, чи будуть замінені на гени 
клітини пам'яті. 


Удосконалений алгоритм оптимізації НІМО-ЗЕ 


Враховуючи розглянуті вище реалізації алгоритмів, котрі реалізують метод мо- 
делювання штучних імунних систем, пропонується наступний алгоритм, що не потребує 
двійкового кодування дійсних змінних та застосовує ряд операторів, відомих, зокрема, 
з генетичних та меметичних алгоритмів. Метою останніх має стати як інтенсифікація 
адаптивного локального пошуку, так 1 якомога повніше вивчення всієї області пошуку. 

Назвемо запропонований алгоритм НІХО-5Е (Нуфбгій Іттипе Меїуотк Оріїтіза- 
поп аїсогійт улйй Зааїу 5еіесіїоп апа КіБопассі 5еатст). Він складається з наступних 
кроків: 

1 Випадково генерується популяція антитіл. Лічильник поколінь встановлюється 
віс0. 

2 Якщо досягнуто наперед заданий максимум кількості поколінь, здійснюється 
перехід до кроку 11; в іншому випадку - перехід до кроку 3. 

3 За оцінкою пристосованості клітин поточного покоління призначається кіль- 
кість клонів для кожної з клітин поточного покоління - оператор селекції. 
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4 Клонування клітин у кількості, що визначена оператором селекції. 

5 Рекомбінація клонованих особин за допомогою оператору ймовірнісного кро- 
соверу. 

6 Розрахунок динамічної ймовірності мутації для клонів та виконання адаптив- 
ного оператора мутацій. 

7 3 певною ймовірністю виконання для кожного клона оператору одновимірного 
локального пошуку за випадковою координатою методом золотого перетину (Фібоначчі). 

8 Стиснення спільної популяції батьків та клонів за рахунок відповідного опера- 


тору до заданого рівня М, . 
9 Знищення з поточного покоління клітин, вік яких досяг заданого значення 1, ,,. 


Їх місце в основній популяції займають клітини, згенеровані випадково рівномірно в 
області пошуку. 

10 Лічильник поколінь 7 - 1-1. Перехід на крок 2. 

11 Виведення поточного покоління, як множини субоптимальних рішень. 

12 Зупинка алгоритму. 

Розглянемо кожен зі згаданих в алгоритмі операторів докладніше. 

Оператор селекції. Для визначення кількості клонів, що припадає на кожну з осо- 
бин поточного покоління, пропонується використовувати оператор, що застосовує метод 
аналізу ієрархій (МАТ) Сааті з відтинанням, вперше застосований у операторах кло- 
нальної селекції еволюційних оптимізаційних стратегій (11). Автор називає свій опе- 
ратор 1 - зрізом множини клітин поточного покоління, що описується рівнянням 


му З Є рІМЛО о кб 1-2М,Ї (6) 
де Йуп 1, - ко - критична пристосованість клітини; ІД) - область пошуку; 


и Ж (Х Й )) - функція належності клітини множині якісних. 


Функція належності клітини визначається наступним чином: 
- рішення нормуються за значенням афінності на відрізок ол; 
и З М ни 
- розраховується рейтинг кожного з них а, -Іч |Ф(Х 5 «|, де Кк - розмірність 


шкали (за умовчанням Х -9); 
- будується матриця попарних порівнянь, для першого рядку якої застосовується 


вираз Б, з а, / а;, а для решти - Б, -б,, /Б,;; 
- функції належності розраховуються як 11 її (Х У )) се зи би 
- КЛІТИНИ, З хз п 


тіп ? 


клонуються у кількості, пропорційній функції належ- 


ності п. -М.: 11 Ж 94 і )) , решта батьківського покоління відкидається. 


Га 
Для ефективної роботи подальших операторів кросовера та мутації, бажано, 
аби особини з високою афінністю продукували достатньо велику кількість клонів. Тому 


рекомендується обирати М, - 10...20-М,,. 


Клональний відбір за МАЇ є однією з ключових відмінностей запропонованого 
алгоритму від решти відомих реалізацій, у яких застосовується турнірний, пропорцій- 
ний («рулетковий») чи рівномірний принципи розподілу кількості клонів між батьківсь- 
кими особинами. Крім того, що даний метод, на відміну від решти, є математично 
обгрунтованим, його практичне застосування забезпечує кращу збіжність при вирішенні 
типових задач у порівнянні з іншими операторами селекції. 

Оператор кросоверу. Іншою ключовою відмінністю запропонованого алгоритму 
є застосування до клонів адаптивного оператору кросоверу, що передбачає обмін ге- 
нетичною інформацією. Оскільки було запропоновано не кодувати дійсну інформацію 


106 «Искусственньгй интеллект» 2013 Мо4 


Алгоритм моделювання штучної імунної системи з селективним оператором... ож 


двійковою, то для рекомбінації генів пропонується застосування модельованого двійко- 
вого кросоверу (5ити/іагеа Біпагу сгоз85оуег, 5ВХ) (1101, що дозволяє імунним клітинам 
виконувати склеювання векторів дійсних координат аналогічно двійковим векторам 
в генетичних алгоритмах. 

Оператор 5ВХ передбачає три параметри, що задаються користувачем: 

- Ру - ЙмОВІрність, що для клона буде виконано оператор кросоверу (зазвичай 


сгозя 
має значення, близьке до 0,5); 

- Р. , - Ймовірність того, що по даній координаті виконуватиметься схрещування 
(може варіюватися в широких межах від 0,4 до І - в останньому випадку кросовер вико- 
нується по всіх координатах одночасно); 

-- індекс варіації п - показує, як сильно нащадок схрещування має бути схожий 
на окремого батька (значення 0 та І з великою ймовірністю дадуть нащадка, схожого 
на одного з батьків, значення 0,5 - максимальне змішування генів). 

Якщо випадкове число, рівномірно розподілене між 0 та 1 (тут 1 надалі позначає- 


мо (0,1)) менше ймовірності Р. ,, для двох випадкових клонів з номерами ) та (, 


СРО55 


а інше випадкове число (/ (031) менше Р 


іпа ? 


то поточна координата і клонів піддається 
схрещуванню за формулою 


РИН: 
хао5 вн вина - вх), (7) 


; 1 
Ха З 5-й т Вже. т (І цій Зеео 
де й - так званий ступінь схрещування, що обчислюється за формулою 
(с «(ол ))пт, у ц(ол)« 4; 
а. (3) 
вче и 
1 їзп 1 
нн хору ЗА , й С(01)»-- 
зе) у (01) а 
Параметр а визначається через ступінь відмінності )-ї та К-ї клітин по ко- 


ординаті і, асаме у -2 


пи Рану - тах Р Хр ) пи (о Хі ї Х піп, і 5 
тах (х пра )- тіп ЕЕ 
Оператор мутації. Головним оператором, що відповідає за пошук рішення, 
в запропонованому алгоритмі, як і в більшості реалізацій методу штучних імунних 
систем, залишається оператор мутації. Пропонується адаптивний оператор, який ви- 
падковим чином застосовує наступну мутацію 
Ге 5 . 
Хі е Хі -б; збо хі ге П :п|, (9) 
де випадкова складова 6, може визначатися на основі Гаусового розподілу згід- 
но з (2) заформулою 


5, - 0,1-М(01), (10) 
або поліноміального розподілу за формулою 
бок бо-окаунря зі, й ЖЕ , ат) 


бо з 


і 


1 
1-6-(-п)н20,-05мигУєто1, оті 


2 
де г, - С (031) - випадкове число, рівномірно розподілене для кожної з коорди- 
нат й о і) - відносне положення поточної клітини в області 
Х 


тах, і Х тіп, і 
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пошуку по і-й координаті; ш - показник щільності поліноміального розподілу (більші 
значення відповідають меншому розсіву мутації). 

Отримуючи на вході традиційні параметри Р" - мінімальний рівень мутації, 
тиї мае -- відносне значення рівня мутації (частини генів, які попадуть під зміну опера- 
тором) та Т.,,, - максимальну припустиму кількість поколінь, алгоритм на кожній іте- 


рації розраховує ступінь мутації генів 


тіп ЗІ й т 
дай на од тика с, Дб» а ее 
тиї (Іеуеі - па М 
тіп З Т 
у зва 2 її ра-те- 
2 


де г - номер поточного покоління (птерації). 

Якщо при переборі координат випадкове число менше ти! (Іеуе! , координата 
піддається тій чи іншій мутації, інакше - залишається без змін. 

Рекомендується задавати Р'"" -- на рівні 0,2...0,5,а тиг таїе - нарівні 0,5...1. 


тиї 
Тоді незалежно від обраного часу роботи алгоритму на початкових ітераціях мутація 
буде виконуватися для всіх клонів, на кінцевих - менше, ніж для половини. 

Вибір - яку саме мутацію виконувати: широку за (10) чи стислу за (11) - здій- 
снюється на основі синтетичної оцінки поточної клітини в масштабі оцінок корисності 
на поточній ітерації по всій популяції. Остання визначається для кожної клітині інди- 
відуально за формулою 


Ї Фах їй Ф; 


Тиах / Фах " Фиіп 
де у - 0,8 - 0,95 - максимальна припустима доля поліноміальних мутацій в за- 
гальній кількості. 

Якщо при зверненні до клона випадкове рівномірне число виявиться меншим 
ти! 5еЇесі , то для генів клона виконуватиметься лише мутація (10), інакше - (11). 

З аналізу (13) зрозуміло, що ширша мутація за Гаусом застосовуватиметься для клі- 
тин з низькою афінністю й частіше -- на початкових кроках оптимізації. Кількість мутацій 
за Гаусом зменшуватиметься лінійно з кількістю ітерацій та пропорційно афінності клі- 
тини на користь поліноміальної мутації. Слід звернути увагу, що вибір типу мутації здій- 
снюється для всієї клітини, а лише потім перевіряється, чи буде виконуватись ця мутація 
для кожної з координат, за (12). 

Оператор локального пошуку. Аби інтенсифікувати локальний пошук рішення 
і зробити його ще більш цілеспрямованим, до клона, який пройшов кросовер 1 мутацію 
з ймовірністю Р, - 0,1...0,25, застосовується оператор одновимірного локального пошуку 


1 
Пиво бен 1 , 7-ДП1:М.), (13) 


методом золотого перетину (Фібоначчі), який складається з наступних операцій |121: 
- рівноймовірно обирається одна з координат К - П 5 п| Е 
- рівноймовірно обирається напрямок пошуку -- до нижньої границі області по- 
і С; ре С. 
шуку тіп Х; чи верхньої тах Х |; 
- здійснюється одновимірний локальний пошук у напрямку від поточного роз- 
; С; б С 
ташування клона у просторі Х;" до границі області пошуку. 
Якщо пошук покращує афінність рішення, знайдене рішення зберігають. 
Оператор стиснення. По закінченні усіх пошукових операторів клітини нового 
покоління оцінюються за якістю. Разом з попереднім поколінням вони утворюють конф- 
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ліктну множину, з якої, використовуючи принцип (1), відбирають У, клітин з найкращою 


афінністю, що лежать не ближче, ніж у радіусі "г одна від одної. Для цього виконуються 
наступні операції: 
- сортування клонів за зменшенням афінності; 
- циклічно від першої клітини до спустошення множини ще не розглянутих; 
- вилучення усіх клітин, близьких до поточної і гірших за неї. 
Слід наголосити також на необхідності вилучення після стиснення клітин, що мають 
вік, більший за визначену межу 7, Пропонується знищувати клітину З 1, х - 10.20 іте- 


рацій, незалежно від її афінності. Подібні клітини замінюються новими, у яких кож- 
на клітина генерується за (3). 


Аналіз результатів роботи алгоритму 


Для порівняння роботи запропонованого алгоритму з відомими реалізаціями 
методу моделювання штучних імунних систем та з іншими евристичними методами, 
що традиційно використовуються для вирішення задач багатовимірної оптимізації, було 
проведено тестування на ряді стандартних багатоекстремальних функцій. Також було 
порівняно роботу запропонованого алгоритму із відомими раніше при вирішенні практи- 
чної задачі - навчання нейронної мережі прямого розповсюдження. 

На рис. 2а представлений знімок екрана результатів виконання алгоритму при 
мінімізації функції Растрігіна Е(Х) «20-х ях; - 10 сов(длх,) - 10 со5(2лх,). Для 
показовості на графік виведено не саму цільову функцію, а зворотну їй величину -- 
афінність клітин. З рис. 2 легко можна побачити основну особливість методу моде- 
лювання імунних систем - отримання в результаті не одного рішення, а одразу пев- 
ної множини, до якої входять як глобальний оптимум, так і найближчі за афінністю 
локальні оптимуми. 


Время виполнения 0.9 сек. (Допустимоє время вюполнения 2 30 сек.) 
Вьполнено 143 итераций  (Допустимое число итераций - 10000 ) 
Количество вичислений целевой функции з 10024 ло 


Ениге 1 
РОРОТГАТІОМ СЕІІ5 ЕІТМЕЗ5 - | Біе вії Мем Іпеегі Тосів Фезкіор Міпфом Неїр 


ря пбень|звіч«09ш4-1а|0 


-0.000303681. 
-0.994960006643824 
-0.994961682318685 


-0.994965004189075 
-0.994972809878648 


-1.990183763175340 
-3.979835314221820 


300 і 350 що 
6) 


Рисунок 2 - Тестування запропонованого алгоритму на функції Растрігіна: 
(а) результат мінімізації п - 2; (б) зміна афінності клітин в часі п - 20 


На рис. 26 показані залежності афінності кращої клітини (безперервна лінія) в по- 
колінні (безперервна лінія) та середньої афінності популяції (пунктир) від номера іте- 
рації при мінімізації функції Растрігіна 20 порядку. Графік, звичний для евристичних 
алгоритмів, має суттєву особливість: навіть найкращі рішення, якщо не були змінені 
протягом певного терміну тривалості життя, відкидаються й замінюються випадковими. 
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У табл. І наведені порівняльні результати роботи алгоритмів при розв'язанні 
типових задач. У якості задачі навчання нейронної мережі була використана задача іден- 
тифікації механічних властивостей сортового прокату від його хімічного складу у скоро- 
ченому та повному варіанті. 


Таблиця 1 - Порівняння роботи алгоритмів при розв'язанні 
р р р рир 


Тестова Розмірність Метод розв'язання задачі 
функція задачі ВС-СА НІА НІХО-8К 
(задача) (кількість Показники якості реалізації алгоритму 
змінних) Б. га я й га й г є я 
Растрігіна 20 0 3,472 | 44210? 0 6,626 | 2,410? 0 1,929 | 105107 
100 1,999 | 15,94 | 1,62107 | 0,999 | 13,33 | 8,610? 0 2,616 | 2,0х10? 
Еклі 20 0 0,142 | 2, 0210? 0 1,086 | 72210? 0 0,007 | 81510" 
100 1,258 | 21,42 | 2,12107 | 1,314 | 11,4 | 942107 | 1,101 | 7,240 | 21210? 
Навчання 72 4,974 | 16,37 | 142109 | 2.994 | 11.33 | 242107 | 2,039 | 10,20 | 2,1-107 
нейромережі 144 14,94 | 39,52 | 51810" | 13,99 | 21.07 | 142107 | 12,07 | 16,75 | 41510? 


Ми є ж зо 5 і -. ха 
В табл. 1 прийняті позначення: / -- мінімальне значення цільової функції, знайде- 
не методом за 10 повторних запусків; /  - статистичне значення оцінки математичного 


очікування мінімального значення функції при випадковому запуску; 7 - статистичне 
значення оцінки математичного очікування кількості викликів цільової функції до за- 
вершення алгоритму при випадковому запуску. 

Як видно, запропонований алгоритм знаходить глобально оптимальні рішення 
частіше від аналогів, при цьому витрачаючи значно менше часу. На особливо склад- 
них задачах, як, наприклад, навчання нейронної мережі на реальних даних, виграш часу 
може становити майже два порядки. 

Аби зрозуміти ефективність кожного з операторів, що складають алгоритм, зве- 
демо у табл. 2 частоту, з якою вони покращують значення афінності клітин. 


Таблиця 2 - Ефективність операторів запропонованого алгоритму 


Оператор Мутація Поліноміальна му- Локальний по- Кросовер 
Гауса тація шук 
Випалки покращення дитеової 23,339 417990, 17,239 114590 
функції, б 


Слід відзначити, що отримані дані вповні відповідають теоретичним засадам, адже 
саме мутація є головним інструментом локального пошуку в методі моделювання штуч- 
них імунних систем. Відсутність у табл. 2 запропонованого оператора селекції пояснюєть- 
ся тим, що він не призводить до покращення рішень, а лише обирає кандидатів для кло- 
нування. 


Висновки 


Аналіз основних сучасних алгоритмів, що реалізують метод моделювання штучних 
імунних систем показав його переваги перед іншими методами багатовимірної умов- 
ної оптимізації у дійсному просторі. Запропоновано гібридний адаптивний імунний 
алгоритм, що використовує оператор клональної селекції, кросовер, адаптивну мута- 
цію і обмежений локальний пошук. 

Селекція здійснюється шляхом оцінки пристосованості рішень з використанням 
методу ієрархії Сааті. Запропонований оператор селекції ефективніший, ніж традиційні 
оператори, крім цього, він є єдиним математично обтрунтованим. 
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Парний адаптивний кросовер клітин-клонів забезпечує інтенсифікацію локального 
пошуку, дозволяючи поєднувати кращі гени двох вдалих рішень з попередніх поколінь. 
Аби підвищити спрямованість пошуку, до випадкової координати клона застосовується 
оператор одновимірної оптимізації методом золотого перетину (Фібоначчі). 

Оператор покоординатної мутації застосовує нормальний або поліноміальний 
закони розподілу ймовірностей. При цьому, якщо клон породжений кліткою з висо- 
кою афінністю, вищою є ймовірність «вузької» поліноміальної мутації, що відповідає 
інтенсивному локального пошуку. Клон, породжений клітиною зі слабкою афінністю, 
має більшу ймовірність «широкої» мутації з відхиленням за законом Гауса. 

Результати дослідження показують високу ефективність запропонованого алгорит- 
му для стандартних цільових функцій, що мають розмірність до 100, а також при вирі- 
шенні завдань навчання нейронних мереж прямого поширення. 

Основні переваги запропонованого алгоритму полягають у тому, що він залишаєть- 
ся ефективним при зростанні розмірності задачі, знаходить не одне рішення, а їх множи- 
ну (альтернативи) та використовує значно менше часу (на порядок) для порівняного 
розв'язання задачі. 
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Тие Аїєогійт ої Атійсіаї Іттипе Зузіет 5ітиатоп улй 5ааїу 


Зеіеспоп Орегаїог апа Опе-аїтепбіопа! Госаї 5еагсй 

Іл Фе агіїсіе ре Базіс тодегп аїсогійрта Бас пріетепі ре агіїйста! шитипе 5узіетя 
зіпиіабоп пефод 15 сопзідегей. А І8о адуапіаєєз апа Ффзадуапіаєез ої (різ аїбогіїитя аге апа- 
Іугед. 

Озіпє а рубгій адаріїує іпатипе аїсогійрга Фаї іпсТтодез ре сіопа! 5еЇесіїоп, сгоз5оуег, 
планоп апа Штікед аЧарбуе ІЇоса! 5еагсб 15 ргорозей аз а пему аЇ ог. 

Зеіеспоп 15 сатіед оці бу ез'штайпо бе 50Їиоп8 заткабійсу цзіпє бе Заагу Апаїуйс Ні- 
егагспу Ргосе85 (АНР). ТРе Бієрег 15 Пе айбпісу ої ре сеї!!5, мпісі. Фебсгібеє ре 50Їшцйоп, Пе 
тоге питкег ої сіопез, Ше сеї! млі! ргоутйе їБе пехі гепеганоп. 

Тре іпуезцєаноп 5ром/8 Шаї Фе ргорозей 8еіесйоп орегагог 15 плоге еНесіує ап (Ре іга- 
дїіопа! орегаїог8, 5исі а5 ре їоцгпатепі, гоціейе апа шпібогта ур сис-ої. Кигфегтоге, із 
зеіеспоп теїрод 15 опіу плапетайсаПу уц5ійНеай. 

Туміп адарбуе сеї! сіопе8 сгоз5омег ргоміде5 ап іпіепзійсаноп ої Іосаї зеагсі апа аПоху8 
їо согабіпе Фе єепе8 ої бо Ше Без! зиссе58їйиі десізіоп5 ої ргеміоця сепегацопя. 

То аргоуе Ше Юосия8 ої'їБе 5еагсп Ше опе-дайтепяіопа! орійтігайоп орегаїог базед оп Фе 
теїодй ої соідеп гапо (патеа Еібопассі) 15 аррпед їо а гапдота соогаїпаїте ої сіопе сеїіз. 

Тре сгоз5омег апа Іосаї 5еагсі орегаїогя аге арріїед їо Ше сіопе сеї8 ул а сегіаїп ргоб- 
абШу, уурбіс сап а5о спапее аЧарбуєїу дигіпе Бе могк об'Фе аїсогійт. 

Соогаїпаїемуіве пликайоп орегаїог уутї ре погта! апа роїупотіа! ргобабійсу дїбітіби- 
бопз аррііез. Могеоуег, 18 а сіопе 15 бепегаїеа Бу а се! уутії Бір аїйпігу, патоуу роіупотіа! 
плианоп 15 плоге ПКеїу ог її. Її соптезропав (о іпіепятує Поса! 5еагср. Іа сіопе 15 сепегагед Бу а 
сеї млі а муеак аїбпігу, а Бгоад пашиіаног улії ре Сацязіап дзілїбибоп 15 плоге (гедїБіе. 

Тре гевийз 5пому Фе рієб еййсіепсу ої'їБе ргорозеа аїсогійрт бог 5гапдага обуесіує Нлпс- 
йоп5 уліб Іе55 Шрап 100 дйтепзіопе, а8 укеі а8 бог 50Іміпє Ше ргобіегт ої Меедїогуага пеига! 
пебмогкя ітаїпіпе. 


Стаття надійшла до редакції 27.06.2013. 
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